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Abstract 

We present a new calculation of the B — > tt form- factor / + , relevant for the measurement 
of \V u b\ from semileptonic B — > 7r transitions, from QCD sum rules on the light-cone. The 
new element is the calculation of radiative corrections to next-to-leading twist-3 accuracy. 
We find that these contributions are factorizable at 0(a s ), which lends additional support 
to the method of QCD sum rules on the light-cone. We obtain /+(0) = 0.26 ± 0.06 ± 0.05, 
where the first error accounts for the uncertainty in the input-parameters and the second 
is a guesstimate of the systematic uncertainty induced by the approximations inherent in 
the method. We also obtain a simple parametrization of the form-factor which is valid 
in the entire kinematical range of semileptonic decays and consistent with vector-meson 
dominance at large momentum-transfer. 



submitted to JHEP 



*E-mail: Patricia. Ball@cern.ch, Patricia.Ball@durham.ac.uk 
^E-mail: zwicky@physik.unizh.ch, Roman. Zwicky@cern.ch 



1. The experimental programme of the dedicated B-factories BaBar and Belle will con- 
tribute to unravel structure and size of flavour- and CP-violation. The key-observable 
is the unitarity-triangle of the CKM-matrix whose overdetermination will help to answer 
the question whether there are additional sources of CP-violation not present in the SM. 
Overdetermination means independent measurements of sides and angles from different 
processes. One of the sides is determined by the CKM-matrix element \V u b\, whose pre- 
cise measurement is certainly a challenging task due to the smallness of the corresponding 
branching ratios. The method of choice is to measure it from semileptonic tree-level decays 
b — > uev where contamination by new-physics effects is expected to be small. The main 
complication in this measurement is, as a matter of course, QCD effects whose calcula- 
tion from first principles is highly challenging. Inclusive semileptonic decays are usually 
treated in heavy quark expansion and become the more delicate the more accurately ex- 
perimentally necessary cuts are taken into account (which necessitates the inclusion of 
other potentially large scales and calls for threshold- and soft-gluon resummation, cf. |]). 
The alternative is to study exclusive decays where, owing to the nonrenormalization of 
vector and axialvector currents, both perturbative and nonperturbative QCD effects are 
neatly encoded in form-factors, depending on only one variable, the momentum-transfer 
to the leptons. The simplest such decay-process, involving only a single form-factoiQ, is 
B — > 7r£b>£, which hence has received fair attention in the literature. First experimental 
results are available from CLEO 0. The most precise calculation of the form-factor will 
sans doute finally come from lattice-simulations; presently, however, the main attention 
of lattice-practitioners appears to be directed not so much to obtaining phenomenologi- 
cally relevant results, but rather to controlling lattice artifacts and approximations like 
discretization errors and the explicit breaking of chiral symmetry with Wilson-fermions, 
particularly relevant for calculating processes that involve pions; another problem is how 
to simulate relativistic b-quarks, and how to access the region of phase-space where the 
pion has large momentum. All these problems are presently under intense debate, and we 
refer to Ref. || for reviews and recent research papers. 

Another, technically much simpler, but also less rigorous approach is provided by QCD 
sum rules on the light-cone (LCSRs) [|], ||. The key-idea is to consider a correlation- 
function of the weak current and a current with the quantum-numbers of the B-meson, 
sandwiched between the vacuum and a pion. For large (negative) virtualities of these 
currents, the correlation-function is, in coordinate-space, dominated by distances close to 
the light-cone and can be discussed in the framework of light-cone expansion. In contrast 
to the short-distance expansion employed by conventional QCD sum rules a la SVZ ||, 
where nonperturbative effects are encoded in vacuum expectation values of local operators 
with vacuum quantum numbers, the condensates, LCSRs rely on the factorization of the 
underlying correlation function into genuinely nonperturbative and universal hadron dis- 
tribution amplitudes (DAs) 4> that are convoluted with process-dependent amplitudes T#, 
which are the analogues to the Wilson-coefficients in the short-distance expansion and can 

1 This statement is true only as long as lepton-masses are negligible, i.e. for semi-electronic and -muonic 
decays. 
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be calculated in perturbation theory, schematically 



correlation function ~ (g) (j)^ n \ (1) 

n 

The sum runs over contributions with increasing twist, labelled by n, which are suppressed 
by increasing powers of, roughly speaking, the virtualities of the involved currents. The 
same correlation function can, on the other hand, be written as a dispersion-relation, in 
the virtuality of the current coupling to the B-meson. Equating dispersion-representation 
and the light-cone expansion, and separating the B-meson contribution from that of higher 
one- and multi-particle states, one obtains a relation for the form-factor describing B — > ir. 

The particular strength of LCSRs lies in the fact that they allow inclusion not only 
of hard-gluon exchange contributions, which have been identified, in the seminal papers 
that opened the study of hard exclusive processes in the framework of perturbative QCD 
(pQCD) 0, as being dominant in light-meson form factors, but that they also capture 
the so-called Feynman-mechanism, where the quark created at the weak vertex carries 
nearly all momentum of the meson in the final-state, while all other quarks are soft. This 
mechanism is suppressed by two powers of momentum-transfer in processes with light 
mesons; as shown in || , this suppression is absent in heavy-to-light transitions^] and hence 
any reasonable application of pQCD to B-meson decays should include this mechanism. 
LCSRs also avoid any reference to a "light-cone wave-function of the B-meson", which 
is a necessary ingredient in all extensions of the original pQCD method to heavy- meson 
decays [§, [|, but whose exact definition appears to be problematic fllUfl . A more detailed 



discussion of the rationale of LCSRs and of the more technical aspects of the method is 
beyond the scope of this letter; more information can be found in the literature [IT]. 



LCSRs are available for the B — > tt form factor / + to 0(a s ) accuracy for the leading 



twist-2 contribution and at tree-level for higher-twist (3 and 4) contributions [|12|, [13|, Q . In 



this letter we calculate the leading radiative corrections to the twist-3 contributions. The 



motivation for this calculation is twofold: first, it has been found in [12], [13|, [14]] that the 
tree-level twist-3 corrections are chirally enhanced and sizeable, and amount up to 30% of 
the final result for /+, which indicates that radiative corrections may be phenomenologically 
relevant. Second, the existence of the factorization-formula ([!]) is nontrivial beyond tree- 
level and, to date, it is only for the twist-2 contribution that factorization has been shown 
to hold also after inclusion of radiative corrections. In this letter we show that ([I]), for 
a certain approximation of the DAs (p^ (leading conformal spin), also holds when 0(a s ) 
corrections to the twist-3 contributions are included. 

2. Let us now properly define the relevant quantities. The form-factors / +i o are given by 

(q = pb- p) 



{2 2 ^ 2 2 
777 TH TH TH 

(PB + p)„ - B q2 * qA + B q2 * fo(q 2 ) M (2) 



in semileptonic decays the physical range in q 2 is < q 2 < (mg — m^) 2 . The starting point 

2 For very large quark masses, though, the Feynman-mechanism is suppressed by Sudakov-logarithms, 
which arc, however, not expected to be effective at the b-quark mass. 
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for the calculation of the form-factor /+ in (^) is the correlation function 



i J dV ?s, (7r(p)|T[n7^](y)K^ 75 rf](0)|0) = U + 2p, + ..., (3) 

where the dots stand for structures not relevant for the calculation of / + . As mentioned 
before, for a certain configuration of virtualities, namely ml — p 2 B > 0(AQ CD mb) and 
ntf, — <? 2 > 0(AQ CD m b ), the integral is dominated by light-like distances and accessible to 
an expansion around the light-cone: 

n + (q 2 , P 2 B ) = J2 /W (n Wm)4 n V;^p^iR)- (4) 

As in (|I|), n labels the twist of operators and fim is the (infrared) factorization-scale. The 
restriction on q 2 , m\ — q 2 > 0(Aq CD mb), implies that /+ is not accessible at all momentum- 
transfers; to be be specific, we restrict ourselves to < q 2 < 14 GeV 2 . As Il + is independent 
of /iiR, the above formula implies that the scale-dependence of must be canceled by 
that of the DAs 

In (|4]) we have assumed that n + can be described by collinear factorization, i.e. that 
the only relevant degrees of freedom are the longitudinal momentum fractions u carried by 
the partons in the n, and that transverse momenta can be integrated over. Hard infrared 
(collinear) divergencies occurring in should be absorbable into the DAs, as discussed 
in detail in Ref. [15]. Collinear factorization is trivial at tree-level, where the b-quark mass 
acts effectively as regulator, but can, in principle, be violated by radiative corrections, by 
so-called "soft" divergent terms, which yield divergencies upon integration over u. Such 
terms break for instance factorization in non-leading twist in the treatment of nonleptonic 
B-decays a la BBNS 0; it is thus instructive to see what happens in the simpler case of the 
correlation function ([3]), where the convolution involves only one DA instead of up to three 
in B — ► 7T7T. To anticipate the result: we find that factorization also works at one- loop 
level for twist-3 contributions and that there are no soft divergencies. 

There are two two-quark twist-3 DAs of the it, <j) p and (j) a , which are defined as 

(0\u(x)[x, -x]ij 5 d(-x)\7r~(p)) = /4(a*ir) / du e l & x (f> p (u, fi m ) , (5) 

J 



(0\u(x)[x, -x}a a ^ 5 d(-x)\7T (p)) = ~ (il(pm)(PaXp - ppx a ) I du e <px c/) a (u, /i m ) (6) 

3 Jo 

with /^j.(/x IR ) = fTTm 2 /(m u + rrid) (hir) and £ = 2u — 1; u is the longitudinal momentum- 
fraction of the total momentum of the 7r carried by the (d-) quark. The Wilson-line 

[x, -x] = P exp[2ig ( dtx^^t - l)x)\ 

^ 

ensures gauge-invariance of the nonlocal matrix-elements. Exploiting conformal symmetry 
of massless QCD, which holds to leading-logarithmic accuracy, one can perform a partial- 
wave expansion of the DAs in terms of increasing conformal spin, which amounts to an 
expansion in Gegenbauer-polynomials and fixes the functional dependence on u; the am- 
plitudes of the partial-waves are of nonperturbative origin and can be related to local 
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Figure 1: Some of the diagrams contributing to T^f'^ in one-loop order. The double line 
denotes the b-quark propagator, the single lines denote the light u- and d-quark propaga- 
tors. 7 M and 75 are the weak and the B-vertex, respectively. There are two more self-energy 
and one more vertex-correction diagrams. 



hadronic matrix-elements by virtue of the QCD equations of motion as discussed in detail 
It turns out, in particular, that the two DAs <p a and <p p are not indepen- 
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IJ, 0. 



dent, but mix with each other and the twist-3 three-particle DA T that parametrizes the 



matrix-element (0|w(x)[x, vxja^^gGap^x^vx, —x]d(—x)\ir~(p)). 
when calculating radiative corrections to H+ to twist-3 accuracy, 
three DAs, 



Thus, for consistency, 
one has to include all 



n 



(3) 



T ®T, 



(T) 



H 



and it is only in the sum of these three terms that hard infrared divergencies and the 
scale-dependence are expected to cancel. The analysis of (16[ has shown, however, that the 
two-quark DAs are very well approximated by the lowest partial-wave, i.e. the one with 
smallest conformal spin, and that mixing with T sets in only at higher conformal spin. 
In calculating radiative corrections to 11+ we thus restrict ourselves to leading conformal 
spin, i.e. the so-called asymptotic DAs, and use jlBj 



(f><j(u) = 6u(l - u), 4> p (u) 



r = 0. 



(7) 



3. Some of the diagrams contributing to T# to one-loop order are shown in Fig. |l| The 
light quarks are massless and have momenta up and up = (1 — u)p, respectively. They 
are projected onto the desired DA by closing the trace with an appropriate projection 
operator V, which is just — ^275/4 for <fr p and involves a derivative in p for (j) a . The 
calculation is performed in dimensional regularization for both ultraviolet and infrared 
divergencies. Carefully distinguishing between the two types of divergent terms, we find 
that the ultraviolet divergencies cancel upon renormalization of the bare b-quark mass in 
the tree-level expression, as they should. The infrared divergent terms, on the other hand, 
do not cancel between /x^ bare 0p( CT ) and 7 i jp( cr )' bare ) separately, but only in the sum of both 

contributions. The renormalized T^' a \u) are regular at the endpoints, i.e. for u — > 0, 
u — > 1, which entails the absence of soft divergent terms. 

As discussed below, the LCSR for / + involves the continuum-subtracted Borel-transform 
B sn \,T H of T H . We calculate it by splitting T H into two terms, T H = yP° le + T^ 1S , where 
T^ 1S can be written as a dispersion-relation in p\ and its Borel-transform is obtained by 
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Figure 2: Radiative corrections to twist- 2 and 3 contributions, in units of a s /(37r), for 
representative input-parameters, as functions of q 2 . Dashed lines: (g) p and (g) </> CT , 
respectively. Solid line: sum of dashed lines. Dotted line: (g) c^ 2 **. The total twist-3 
correction (solid line) is much smaller than the twist-2 correction (dotted line). 

applying Eq. ([5]). T^ ole has a (single or double) pole in s = m\ — up 2 B — uq 2 — * 0; the 
Borel-transforms of these terms are a bit more involved, master-formulae are given in the 
appendix. The final expressions for B su ^T B ' a ^ are too bulky to be presented here. A com- 
pact version of the tree- level expression can be found in In Fig. ^ we compare the 
relative size of radiative corrections to twist-2 and 3 contributions. Whereas the absolute 
value of the two twist-3 corrections is, separately, of roughly the same size as that of the 
twist-2 correction, the only relevant quantity, the sum of both twist-3 corrections, is much 
smaller than the twist-2 contribution, which indicates a good convergence of the light-cone 
expansion also at 0(a s ). 

The result of this calculation is the light-cone expansion of Il + , IT^ , and its continuum- 
subtracted Borel-transform, i? su bn^ c . 

4. Let us now derive the LCSR for / + . The correlation function Il + , calculated for 
unphysical p 2 B , can be written as dispersion-relation over its physical cut. Singling out the 
contribution of the B-meson, one has 

n + = f+(q 2 ) — „ B ^ B ' + higher poles and cuts, (8) 

m B Pb 

where Jb is the leptonic decay constant of the B-meson, fs m \ = 7?7.&(-B|6z75cf|0). In the 
framework of LCSRs one does not use (19) as it stands, but performs a Borel-transformation, 

with the Borel-parameter M 2 ; this transformation enhances the ground-state B-meson 
contribution to the dispersion-representation of 11+ and suppresses contributions of higher 
twist to the light-cone expansion of 11+ . The next step is to invoke quark-hadron duality 
to approximate the contributions of hadrons other than the ground-state B-meson by the 
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imaginary part of the light-cone expansion of so that 

BTL™ = ^ m%f B f + (q 2 ) e- m Bl M ' 2 + ^ ± Jjt Imnf (t) exp(-i/M 2 ) (10) 

and 4ubITf = ^ m%f B f + {q 2 ) e-**'* . (11) 

Eq. (|TTD is the LCSR for / + . s is the so-called continuum threshold, which separates the 
ground-state from the continuum contribution. At tree-level, the continuum-subtraction 
in (11) introduces a lower limit of integration, u > (m 2 — q 2 )/(so — q 2 ) = Uq, in (|j), which 
behaves as 1 — Aq C d / mb f° r large m h and thus corresponds to the dynamical configuration 
of the Feynman-mechanism, as it cuts off low momenta of the u-quark created at the weak 
vertex. At 0(a s ), there are also contributions with no cut in the integration over u, which 
thus correspond to hard-gluon exchange contributions. Numerically, these terms turn out 
to be very small, ~ 0(1%) of the total result for /+. As with standard QCD sum rules, the 
use of quark-hadron duality above s and the choice of s itself introduce a certain model- 
dependence (or systematic error) in the final result for the form-factor, which is difficult to 
estimate. In this letter we opt for being rather conservative and add a 20% systematic error 
to the final result for / + . Another hadronic parameter showing up in ([11]), which actually 
allows one to fix the value of Sq, is fs can in principle be measured from the decay B — ► 
£p£, which, due to the expected smallness of its branching ratio, BR~ O(10~ 6 ), has, up to 
now, escaped experimental detection, fs is one of the best-studied observables in lattice- 
simulations with heavy quarks; the current world-average from unquenched calculations 
with two dynamical quarks is fs = (200 ± 30) MeV |L8]. It can also be calculated from 
QCD sum rules: the most recent determinations |19j include 0(a 2 ) corrections and find 
(206 ± 20) MeV and (197 ± 23) MeV, respectively. For consistency, we do not use these 
results, but replace Jb in ( |TTD by its QCD sum rule to 0(a s ) accuracy, including dependence 
on So and M 2 . For the b-quark mass, we use an average over recent determinations of the 
MS mass, rn 'f ) ^Msi m b) ~ (4.22 ± 0.08) GeV ||20| , |21|| , which corresponds to the one-loop 
pole-mass m^iL-pde = (4.60 ± 0.09) GeV; this is different from the value = 4.8 GeV 
we used in our previous paper, the first reference in jn|. With these values we find 



J'b = (192 ± 22) GeV (the error only includes variation with m& and M 2 , at optimized 
s ), in very good agreement with both lattice and QCD sum rules to 0(a 2 ) accuracy. For 
m b = (4.51, 4.60, 4.69) GeV the optimized s are (34.5, 34.0, 33.5) GeV 2 , and the relevant 
range in M 2 is M 2 ss (4.5-8) GeV 2 . 

The infrared factorization-scale is set to /x 2 R = m 2 B — m 2 p2]| ; the dependence of /+ 
on nm is very small, as all numerically sizeable contributions are now available in next- 
to-leading order in QCD, which ensures good cancellation of the scale-dependence. For 



the 7r-DAs we use for the most part the same expressions as in the first reference of [|13 
except for the new 0(a s ) corrections to Tjj , where we use the DAs given in (^), and for 
the twist-2 DAs, where new analyses of the available experimental data on the 7*77? and 
the 7r electromagnetic form-factor indicate that is closer to its asymptotic form than 
assumed previously and well approximated by p3 | 



(2) (n,/i IR ) = 6«(1 -u) jl + a 2 (/ii R ) (— (2u 



15 ,„ 3 



1 --) 1 (12) 



6 



1 

. 8 
gO. 6 
^+0.4 

0.2 

2 4 6 8 10 12 14 

Figure 3: f+(q 2 ) as function of momentum-transfer q 2 from the LCSR (|TT|). Solid line: 
LCSR for central values of input-parameters and M 2 = 6 GeV 2 . Dashed lines: dependence 
of /+ on variation of input-parameters as specified in the text and 5 GeV 2 < M 2 < 8 GeV 2 . 




with 02(1 GeV) = 0.1 ± 0.1. The two-particle twist-3 DAs (fi p and CT are proportional 
to /x 2 = m 2 / 7r /(m u + rrid), which in the chiral limit equals —2(qq)/f n . Using the "stan- 
dard value" of the quark-condensate, (qq)(lGeV) = (—0.24 GeV) 3 , one has /i 2 (/i 2 R ) = 
0.25 GeV 2 . With the central value for the flavour-averaged light-quark mass, m u d{2 GeV) = 



4.5 MeV, from lattice one has [i 2 (/x 2 R ) = 0.30 GeV 2 , with a quoted uncertainty of 20%. 
In our analysis we use the average /i 2 (/i 2 R ) = (0.27 ± 0.07) GeV 2 . 

With these parameters, we obtain the results shown in Fig. |3|. The form-factor can be 
accurately fitted by 



l-a(q /m B ) + b{q fm B ) 

with /+(0), a and b given in Tab. 1, for different values of m&, so and M 2 . The above 
parametrization reproduces the actual values calculated from the LCSR, for q 2 < 14 GeV 2 , 
to within 2% accuracy. At q 2 = we find /+(0) = 0.26±0.06, when all input-parameters are 
varied within the ranges specified above. This has to be compared to 0.31 ± 0.05 from our 
previous analysis fl3| . At fixed M 2 = 6 GeV 2 , and the central value for nib, rrib = 4.6 GeV, 



we obtain, in this letter, /+(0) = 0.261. Without the new radiative corrections, the result 
becomes 0.294. Using in addition the parametrization of (p^ employed in [13], this becomes 



0.321. And switching to nib = 4.8 GeV and sq = 33.5 GeV 2 , one obtains 0.308, i.e. the 



central value for /+(0) quoted in [|13|]. The tree- level result for the central values of input- 
parameters is 0.247, i.e. the total effect of radiative corrections is below 10%. 

It is instructive to compare the parametrization (|T3"D, obtained at not too large q 2 , 
q 2 < 14 GeV 2 , with the vector-meson pole- dominance approximation valid at large q 2 ~ 
{m B — m^) 2 = 26.4 GeV 2 : here /+ is dominated by the B*-pole, located at q 2 = m 2 B * = 
28.4 GeV 2 , and can be expressed as 

U(q 2 ) = -^- (14) 

1 — z^r~ 



D" 



The residue of the pole, c, can be related to physical couplings as c = f b* 9bb*-k / (2m ) , 
where Jb* is the leptonic decay constant of the B* and Qbb*-k is the coupling of the B* to 
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m b [GeV] 


4.69 


4.60 


4.51 


s [GeV 2 ] 


33.5 


34.0 


34.5 


a 2 (lGeV) 





0.1 


0.2 


/4(/4i)[GeV 2 ] 


0.34 


0.27 


0.20 


M [Gey J 


8 


6 


5 


/+(o) 


0.206 


0.261 


0.323 


a 


2.34 


2.03 


1.76 


6 


1.77 


1.29 


0.87 


g 2 [GeV 2 ] 


14.3 


15.7 


18.5 


Cfit 


0.384 


0.439 


0.523 


}b*9bb*tt 
clcsr — 

2m r* 


0.396 


0.414 


0.430 



Table 1: Fit-parameters of Eqs. and (|T4|), including all variations of input-parameters 
and M 2 . 



the B7r-pair. c can be calculated from LCSRs itself: it is known up to twist-4 at tree-level; 
0(a s ) radiative corrections are known for the twist-2 contribution |2^, |24|. 

We can now try to match the parametrizations (|13|) and ([14]) . To this end, we treat c 
as fit-parameter and require that the transition between both parametrizations be smooth, 
i.e. that at q$ to be fitted, the values of both f+(q$) and its derivative are equal for both 
parametrizations. The resulting values of q$ and c^ t are tabulated in Tab. 1, and the 

< ql and (|14|) for q 2 > 



corresponding form-factors, obtained from plotting for q 2 
are shown in Fig. |j. The last row in Tab. 1 gives the values of clcsr obtained directly from 
the LCSR calculated in 



22, 



IM. The agreement between the direct and the fitted values is 
remarkably good, in view of the fact that the LCSR for c is less accurate than ours for / + , 
as it does not include 0(a s ) corrections at twist-3. Also the values of g 2 are well within 
expectation: sufficiently below the pole on the one hand, but not too small on the other 
hand. Motivated by these results, we suggest a new parametrization of / + in terms of 5 
parameters: Eq. (|13"D for q 2 < q^ and Eq. (0) for q 2 > q 2 ., with the set of parameters 
given in Tab. 1 which comprises the full dependence of f+(q 2 ) on input-parameters and 
the Borel-parameter. 

5. In this letter we have investigated the light-cone sum rule for the form-factor f?^, 
including the calculation of 0(a s ) radiative corrections to the next-to-leading twist-3 con- 
tribution. The calculation has demonstrated the validity of the factorization formula (f|) 
and the absence of soft divergent terms to the considered accuracy, i.e. twist-3 to 0(a s ) 
with the DAs in the approximation of leading conformal spin, Eq. (^). As already found 
in |TJ]], the Feynman-mechanism is the dominant contribution to /+ and hard perturbative 
corrections are numerically small. 

In view of the systematic uncertainties inherent in LCSRs, further refinement of the 
LCSR for f+^ n by including even higher twist contributions or more perturbative correc- 
tions is not likely to increase the overall accuracy. Improvement of the central value of the 
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Figure 4: f+(q 2 ) as function of q 2 in the entire physical range in B — > Treu, from Eqs. (13) 
and fll4"|), including all variations of input-parameters and M 2 . Solid line obtained from 
3rd column in Tab. 1, lower dashed curve from 2nd column and upper dashed curve from 
4th column. 



result may, however, come from a reduced uncertainty of input-parameters: whereas the 
dependence of f+(q 2 ) on is rather small, a reduction of uncertainty in the non- leading 
conformal spin contributions to <fi^ would be useful. As for the approximation of leading 
conformal spin DAs in twist-3, Eq. (0), at tree- level it yields about 90% of the actual result, 
which strengthens confidence that this approximation works similarly well also at 0(a s ). 
The change in central values of f+(q 2 ) as compared to our previous results, Refs. ||13|, [L4], 
is partially due to the 0(a s ) corrections to twist-3 we have calculated in this letter and 
partially due to updated input-parameters. 

The LCSR is valid only for large energies of the tt, i.e. not too large values of q 2 . In 
this letter we have fixed, somewhat arbitrarily, the maximum allowed value of q 2 at g^ ax 
at 14GeV 2 and have parametrized the form-factor by Eq. (|13|). On the other hand, for 
q 2 close to the kinematical maximum allowed in semileptonic decays, q 2 = 26.4 GeV 2 , the 
form-factor is dominated by the close- by B*-pole and can be parametrized by Eq. ([14]). 
The residue of that pole can also be calculated from LCSRs. We have tried to match 
both parametrizations, Eqs. flT3| ) and fll4"|), by requiring smoothness at the matching-point 
q 2 = q 2 , which is a parameter of the fit itself. The resulting values of q 2 are within 
expectations, and the fitted values of the residue agree well with the direct calculation 
from LCSRs, as demonstrated in Tab. 1. This result lends additional confidence to our 
final parametrization of the form-factor, i.e. the combination of Eqs. ( [IB] ) and (TT4]), with a 
total of 5 parameters, which is valid in the complete range of kinematically allowed q 2 in 
B -> vrez/, < q 2 < (m B - m n ) 2 = 26.4 GeV 2 . 



Appendix 

In this appendix we collect formulas for non-standard Borel-transformations. Generally, 
the Borel-transform Bf(P 2 ) of a function f{P 2 ) of the Euclidean momentum P is defined 
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as 

P 2 /N=M 2 fixed 

By B sn i,f(P 2 ) we denote the Borel-transform including continuum-subtraction above the 
threshold So, i.e. if / has the dispersion-representation (here p is Minkowskian) 



f(p 2 ) = fdt 

we define 



P(t) 
t-p 2 



4ub/(p 2 ) = ^ Jdtp(t)e^ M2 . (A.l) 



We need in particular the following transforms (s = m 2 — up 2 — uq 2 ^ 



m — ut — 



(t-^sP r(a + /3) uP(M 2 ) a +p 1 1 v «m 2 y 

from which the Borel-transforms of expressions with additional logarithms are obtained as 
derivatives, e.g. i?s _/3 lns = -^5s"' 3 . Including continuum subtraction, we find 

2 -2 

m — uq / i-22 

n 1 _ e~^~ i / r(i-/?, , , . . 
Bs » h ^-J^pyW)[ n^P) — i° (M - Mo) 

with Uq = (m 2 — q 2 )/(s — q 2 ). For integer j3, the second terms becomes a sum over 
8{u — Uq) and its derivatives. We also give the spectral function p of the general expression 
(m 2 — p 2 )~ a s~ 13 , from which the Borel-transform -B su b including continuum subtraction can 
be obtained using ( |A.1| ): 

/x ^ / ~ 9 x ^ / 9 , sin o;7r 1 
pit) = Q(m 2 - t)0(t - m 2 ^ 



ii (t-m 2 ) a uP(m 2 - if 



+ e(t-m 2 )* in(a + /3)7r 1 



ii (t-m 2 ) a uP{t-fh?Y 
with m 2 = (m 2 — uq 2 )/u. 
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